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ABSTRACT 

Context. Sausage modes are important in coronal seismology. Spatially damped propagating sausage waves were recently observed 
in the solar atmosphere. 

Aims. We examine how wave leakage influences the spatial damping of sausage waves propagating along coronal structures modeled 
by a cylindrical density enhancement embedded in a uniform magnetic field. 

Methods. Working in the framework of cold magnetohydrodynamics, we solve the dispersion relation (DR) governing sausage waves 
for complex-valued longitudinal wavenumber k at given real angular frequencies a>. For validation purposes, we also provide analytical 
approximations to the DR in the low-frequency limit and in the vicinity of a> c , the critical angular frequency separating trapped from 
leaky waves. 

Results. In contrast to the standing case, propagating sausage waves are allowed for a> much lower than oj c . However, while able 
to direct their energy upwards, these low-frequency waves are subject to substantial spatial attenuation. The spatial damping length 
shows little dependence on the density contrast between the cylinder and its surroundings, and depends only weakly on frequency. 
This spatial damping length is of the order of the cylinder radius for cj < 1.5i)Ai /a, where a and UAi are the cylinder radius and the 
Alfven speed in the cylinder, respectively. 

Conclusions. If a coronal cylinder is perturbed by symmetric boundary drivers (e.g., granular motions) with a broadband spectrum, 
wave leakage efficiently filters out the low-frequency components. 
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1. Introduction 


Considerable progress has been made in coronal seis¬ 
mology thanks to the abundantly identified waves and 
oscillations in the structured sol ar atmosp here (for a re¬ 
cent review, see [D^Mportel_&^akmako^_ 2012, and 

Nakariakov & Erdelvil 120091 


also Ballester et al.l 2007, 


Erdelvi & Goosseria 1201 ll for three recent topical issues). 
Equally important is a detailed theoretical understanding of the 
collective wav e modes supported by magnetized cylinders (e.g., 
Robert si 1200 0). While kink waves (with azimuthal wavenumber 
m — 1) have attracted much attention since their measurement s 
with TRACE dNakariakov et al.lll999tlAschwanden et al.|[T999l) . 
sausage waves prove important in interpreting second-scale 
quasi-periodi c pulsati ons (QPPs) in the lightcurves of solar 
flares (see [Nakariakov & Melnikov! 120091 for a recent review). 
Their importance is strengthened given their recent detec¬ 


tion in both the chromosphere ([Morton et al 


photosphere dDorotovic et al .1120141 : [Grant et al .1120151) . 


20121) and the 


Standing sausage modes are well understood. For in¬ 
stance, two distinct regimes are known to exist, depending 
on the longitudinal wavenumber k dNakariakov & Verwichtei 
120051 ). The trapped regime results when k exceeds some 
critical value k c , whereas the leaky regime a rises when the 
oppos i te is true. Both eigen-mode analyses dKonvlova et al.l 
2007 : IVasheshani Farahani et al.l 12014 ) and numerical sim- 


ulations from an in itial-value-problem perspective (e.g., 
iNakariakov et al.l 12012) indicate that the period P of sausage 
modes increases smoothly with decreasing k until reaching 
some Pq in the thin-cylinder limit (ka —> 0 with a being 
the cylinder radius). Likewise, being identically infinite in 
the trapped regime, the attenuation time r decreases with 
decreasing k before saturating at to when ka — 0. Furthermore, 
P 0 is determined prim arily by a/vM , where i: A i is the Alfven 
speed in the cylinder dRoberts et al.l 119841) , with the detailed 
t ransverse density distri butio n playing an im portant role 
dNakariakov et al. 2012 . also IChen et alJ l2015l> . This is why 
second-scale QPPs are attributed to standing sausage modes, 
since a/vAi is of the order of seconds for typical coronal 
structures. On the other hand, the ratio tq/Pq is basica lly 
proportional to the density contrast (e.g.. lKopvlova et alj|2007l) . 
meaning that high-quality sausage modes are associated with 
coronal structures with densities considerably exceeding their 
surroundings. 

Interestingly, the dispersive behavior of trapped modes is 
important also in understanding impulsively generated sausage 
waves (Roberts et alill984 ). When measured at a distance from 
the source, the signals from these propagating waves possess 
three phases: periodic, quasi-periodic, and decay. The frequency 
dependence of the longitudinal group speed u gr is crucial in 
this context. In particular, whether the quasi-periodic phase 
exists depends on the existence of a local minimum in u gr , 
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which in turn depends on the density profile transverse to th e 
structure (lEdwin & Robertsll988tlNakariakov & Roberts! 19951) . 
This analytical ex pectation, extensively reproduced in numerical 
simulations (e.g.j Muraws ki & Robertslll993llSelwa et al .1120041 
iNakariakov et al.ll2004l) . well explains the time signatures of the 
wave tra ins discovered with the Solar Eclipse Coronal Imag ing 
System dWilliams et al.ll20()H 120021: iKatsivan nis et alJ 2003b as 
well as those measured with SDO/AIA dYuan et alJl2013l) . 


We intend to examine the spatial damping of leaky 
sausage waves propagating along coronal cylinders in re¬ 
sponse t o photospheric motions due to, say, granular con¬ 
vection dBerghmans et al. Il996l) . One reason for conduct¬ 
ing this study is that, besides the observations showing 
that pr opagati n g sausag e waves abound in the chromo¬ 
sphere (IMorton et alJ 12012 1. a recent study clearly demon¬ 
strates the spatial damping of sausage waves propagating from 
the photosphere to the transition region in a pore iGrant et al.l 
l2015lt. Another motivation is connected to the intensive int er¬ 

est dTerradas et alJl2010l : lHood et al.l2013HPascoe et al.l2013l) in 
employing resonant absorption to understand the spatial damp¬ 
ing of propagating kink waves measured with the ^Coronal 
Multi - Channel Polarimeter (CoMP) instrument dTomczvk et al.l 
120071 : iTomczvk & Mclntoshl 12009 1. A leading mechanism 
for int erpreting the temporal damping of standing kink 
modes dRuderman & Robertsl 120021 lAschwanden et all 120031 
and references therein), resonant absorption is found to atten¬ 
uate propagating kink waves with a spat ial length inversely pro¬ 
portional to wave frequency dTerradas et al.ll2010l) . If attributing 
the generation of these kink waves to broadband photospheric 
perturbations, one expects that resonant absorption essentially 
filters out the high-frequency components. One then naturally 
asks: What role does wave leakage play in attenuating propagat¬ 
ing sausage waves? And what will be the frequency dependence 
of the associated damping length? 

This manuscript is structured as follows. We present the nec¬ 
essary equations, the dispersion relation (DR) in particular, in 
Sect|2l and then present our numerical solutions to the DR in 
Sect[2 where we also derive a couple of analytical approxima¬ 
tions to the DR for validation purposes. Finally, a summary is 
given in Sect]?] 


2. Problem Formulation 

We consider sausage waves propagating in a structured corona 
modeled by a plasma cylinder with radius a embedded in a 
uniform magnetic field B = Bz., where a cylindrical coordi¬ 
nate system ( r,6,z ) is adopted. The cylinder is directed along 
the ^-direction. A piece-wise constant (top-hat) density profile 
is adopted, with the densities inside and external to the cylinder 
being p\ and p t , respectively (p e < p,). The Alfven speeds, (j Al 
and UAe. follow from the definition va = \jB 2 /4nf>. Appropriate 
for the solar corona, zero-/?, ideal MHD equations are adopted. 
In such a case, sausage waves do not perturb the z-component of 
the plasma velocity. Let 6v r denote the radial velocity perturba¬ 
tion, and 6b r , 6b z denote the radial and longitudinal components 
of the perturbed magnetic field 6b, respectively. The perturbed 
magnetic pressure, or equivalently total pressure in the zero-/? 
case, is then 6p tot = Bdb-JAn. 

Let us Fourier-decompose any perturbed value 6f(r, z; t) as 

6f(r, z; t ) = Re \f(r) exp (tot - fe)]} . (1) 
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With the definition 

, or ,, , to 2 ,, , n 7r 

Bi = ~2 - Bo = — ~ k (-T < ar 8AA ar g/L < -), (2) 

V ~Ae 2 2 

linearizing the ideal MHD equations then leads to (e.g. JCallvi 
11985) 

~ r (rp'tot)' + B 2 Ptot = 0, (3) 

where the prime ' — dldr, and this equation is valid both in¬ 
side (denoted by the subscript i) and outside (the subscript e) 
the cylinder. The solutions to Eq. ([3]) are given by 

Ptot,i = AiMpir), p t ot,e = A e H ( 0 l \p e r), (4) 

where J n and //,) 11 are the «-th-oder Bessel and Hankel functions 
of the first kind, respectively (here n = 0). For future reference, 
we also give the expressions for the Fourier amplitudes of some 
relevant perturbations. 


b z 


Anp lol 

B 


v r 



aT kB ~ 
and b r =- v r . 

to 


(5) 


In addition, the energy flux density carried by the sausage waves 
is given by 


F = 


1 6v r 6b r 

—— (6\ x B) x <5b = 6p tot 6v r r --- Bz. 

47T 47T 


(6) 


The dispersion relation (DR) for sausage waves follows from 
the requirements that the Fourier amplitudes of both the total 
pressure perturbation p tot and the radial velocity perturba¬ 
tion u r be continuous at the cylinder boundary r — a. Note that 
more properly, the continuity of the Fourier amplitude of the 
radial Lagrangian displacement f r should be used in place of 
that of v r . Even though in the static case the two requirements 
are equivalent given that v r = -ito£ r , th e equivalence will not 
be prese nt when axial flows exist (e.g..|Goossens et alJl!992c 
iLi et alJl20H . The DR reads (e.g.. lCallvill986ir ~ 


IA _ /i(pjfl) 

Be Jo(Bi a ) H[ l \p e a)' 


Throughout this study, we examine only the solutions that cor¬ 
respond to the lowest critical longitudinal wavenumber in the 
trapped regime. Unless otherwise specified, we solve Eq. 0 by 
assuming that the angular frequency to is real, and the longitudi¬ 
nal wavenumber k is complex (k - k^ + ik\). For these propagat¬ 
ing waves, the energy flux density averaged over a wave period 
2nlto is 

(F r ) = ^Re (pioiK) exp(-2kiz), (8) 

( F z) = ^Re(-^M;jexp(- 2 kiz), ( 9 ) 


where the subscripts r and z denote the radial and longitu¬ 
dinal components, respectively. Furthermore, /* means taking 
the complex conjugate of some complex-valued /. (See Ap¬ 
pendix A for a derivation of Eqs. @ and ©.) 

Before proceeding, we note that a comprehensive study 
was conducted recently by iMoreels et alJ ( 2015 1 hereafter 
M15) to work out the specific expressions for the energy and 
energy flux densities for both fast and slow sausage waves. 
That study differs from ours primarily in the objectives. To 
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enable a calculation of the energy content of sausage waves 
in a variety of solar structures, M15 employed the specific 
expressions for the eigen-functions and focused on trapped 
waves. However, the aim of the present study is to exam¬ 
ine the spatial damping of propagating leaky waves. For this 
purpose, the specific expressions for the energy flux densities 
are not necessary. Rather, we only need to evaluate the sign 
of (F r ) at large distances from the cylinder to offer a phys¬ 
ical explanation for the spatial damping. Likewise, the sign 
of (F z ) is necessary to address whether the spatially damped 
waves can impart their energy upwards. We would like to 
stress that, while the energetics of leaky waves cannot be ad¬ 
dressed with an eigen-value-problem approach (see our dis¬ 
cussion towards the end of Sect. [3), evaluating the signs of 
(F r ) and (F z ) makes physical sense. A similar conclusion was 
reached in M15 for the trapped waves at the critical axial 
wavenumber, below which trapped waves are no longer al¬ 
lowed. Regarding the technical details, in the present study 
the energy flux density, Eq. ©, is equivalent to the Poynting 
flux given the cold MHD approximation. In other words, the 
contribution due to thermal pressure, given by (T) in Eq. (4) 
of M15, vanishes here. 


3. Results 

Figure [T] presents the solutions to the DR for two different den¬ 
sity ratios, one mild (pjp c — 4, the red curves and symbols) 
and the other rather large (f> JfP - 100, blue). The real (the filled 
dots) and imaginary (open) parts of the longitudinal wavenumber 
k are presented as functions of the angular frequency u>. The ver¬ 
tical dash-dotted lines separate the trapped regime (where k\ = 0) 
from the leaky one (k t 4 0), and correspond to op = k c VAe with 
k Q = 2.4048/(a yjpi/p e - 1) JNakariakov & Verwichtel 1200 51 
Let us leave the physical interpretation of the results till later, 
and for now provide a validation study on the numerical solu¬ 
tions. To this end, the dispersion curves in the neighborhood of 
at c are amplified in Fig. EH- Assuming that |Aw| = \<±>- w c | «; at c , 
one can derive an analytical expression. 


1 + " { 

In 

k 2 c a 2 / 

~ V 

A oj _ A k \ 
co c k c ) 

+ 2 r-| 

} k c 

1 + 4 

7T 

{hr 

’ k 2 a 2 

2 

( Aco A k' 
\0Jc kj 

+ 2y — 1 j 

at c 


( 10 ) 


where A k — k - k c . Interestingl y, Eq. (flOt agrees with 
Eq. (16) in lVasheghani Farahani et all (i2014l) even though stand¬ 
ing modes with real k and complex-valued at are examined there. 
The primary difference between the two is the appearance of the 
Euler constant y = 0.577, which we found is necessary to retain 
when the logarithmic term does not substantially exceed unity. 
In Fig. [TJj, the solid curves represent the solutions to Eq. ( I I Ob . It 
can be seen that they well approximate the solutions to the full 
DR given by the dots in both trapped and leaky regimes. 

Another portion of the dispersion curves that is analytically 
tractable is when at —> 0. In this case, the DR can be approxi¬ 
mated by 


p x a — arctanf———) . (11) 

\t*e + m) 

To arrive at Eq. CD. we have assumed that |yUjfl|, |/r e fl| » 1, 
which can be justified a posteri. The real (imaginary) part of 
the solution to Eq. (fill is presented by the curves in Fig. EH 
(Fig. EH), and is found to agree well with the solutions to the 



co a/' t' Ai co a/v M 


Fig. 1. Dependence of the real ( k r, the solid dots) and imaginary (k\, 
open dots) parts of the longitudinal wavenumber on the angular fre¬ 
quency at. Two density ratios pjp e = 4 and 100 are examined, and are 
given by the symbols and curves in red and blue, respectively. The ver¬ 
tical dash-dotted lines correspond to the critical angular frequency at c 
that separates the leaky (to the left of w c ) from trapped (right) regimes. 
Figure [TJi presents an overview of the dispersion curves, while FigEH 
(Figs. EH and EH) examines the portion where at is close to at c (at ap¬ 
proaches zero). The dots are found by numerically solving the full dis¬ 
persion relation (DR, Eq. 0) by assuming a real at but a complex¬ 
valued k = &r + iki, while the curves represent the solutions to the 
approximate DR (see text for details). 



C3 

-y 


Fig. 2. Dispersion curves of (a) propagating and (b) standing sausage 
waves. In (a) the real (k R ) and imaginary (k\) parts of the longitudinal 
wavenumber are given as functions of the real angular frequency a>. In 
(b) the real (m R ) and imaginary (wi) parts of at are shown as functions 
of the real wavenumber k , which is given by the vertical axis though. 
In both panels, the solid (dashed) curves represent the trapped (leaky) 
regime. Two density ratios (pjp t = 4 and 100) are examined as repre¬ 
sented by the red and blue curves, respectively. 


full DR, represented by the dots. It is also clear that in the low- 
frequency portion the solutions to the DR for the two density 
ratios are close to each other, despite that the ratios differ sub¬ 
stantially. This is understandable from Eqs. (ITU and (0, because 
/./? and p 2 both approach -k 2 when at approaches zero, therefore 
having little dependence on the Alfven speeds or densities. 

Figure|2ji presents once again the dispersion curves for prop¬ 
agating sausage waves, where the real and imaginary parts of k 
are given as functions of at. For comparison. Fig. EH presents 
the dispersion curves for standing waves, obtained by solving 
Eq. (0 for complex-valued at (at = wr + iat i) at given real 
values of k. Note that here k is given by the vertical axis, and 
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-u)\ is plotted instead of toi because a>\ < 0. In both cases, we 
examine two density ratios with pjp e being 4 (the red curves) 
and 100 (blue), and present the trapped (leaky) regime with the 
solid (dashed) curves. One can see that in the trapped regime. 
Figs. Eh and [2j> agree exactly with one another, as expected 
for the situation where both a> and k are real. In the trapped 
regime, however, some distinct differences appear. For stand¬ 
ing waves, with decreasing k, the real (imaginary) part of the 
angular frequency wr (oji) decreases (increases in magnitude) 
and saturates when k approaches zero. This_k-dependence of 
both u>r and ut; is wel l understood (e.g., Kopylova et al.l 120071 : 


iNakariakov et al][2012t IVasheghani Farahani et al.ll2QI4l) . It suf¬ 

fices to note that the temporal damping clearly depends on p;/p e : 
the higher the density ratio, the slower the temporal damping. In 
fact, for large pi/p e , a simple but reasonabl y accurate exp r essio n 
exists, namely, |u>i|/w R « (7r/2)(p e /pi) (IKonvlova et al.l 12007 ). 
When examining propagating waves. Fig. [2ji shows that the DR 
permits waves with to much lower than to c , in contrast to the 
standing case. In addition, for a> < v^J a, neither k K nor k\ shows 
a significant dependence on the density ratio. This is particularly 
true for k\, and has been explained in view of Eq. (ITTI) . 

Why should the results be different if one simply changes 
from one perspective, where a DR is solved for complex to as 
a function of real k, to another standpoint where the same DR 
is solved for complex k as a function of real to"! As has been 
discussed in detail by Tagger et al.l ( 1995 . hereafter TFS95), 
indeed the two perspectives have a close relationship when 
wave attenuation is weak or absent. However, some consid¬ 
erable difference between the two may result when strong 
attenuation takes place. An example for this can be found 
in TFS95 where the authors examined Unear Alfven waves 
in a partially ionized gas where ions and neutrals are im¬ 
perfectly coupled. Solving the relevant DR for complex a> 
at real k yields that to may be purely imaginary in certain 
ranges of k, in other words, the waves may be overdamped 
and non-propagating. In contrast, solving the same DR for 
complex k at real to yields that the real part of k is always 
non-zero, meaning that the waves are always propagating. 
This led TFS95 to conclude that how to choose a perspective 
depends on how the waves are excited. One chooses complex 
k and real to to examine the spatial variation of the waves 
excited at a given location with a given real frequency. On 
the other hand, one chooses complex to and real k to follow 
the temporal variation of the waves in response to perturba¬ 
tions initially imposed in a coherent way over many wave¬ 
lengths. These discussions are also valid if one examines the 
differences in Figs. Eh and [2}). In particular, the absence of 
standing waves with frequencies below a certain value (see 
Fig- Eh) can be understood given that in view of Fig. [2h, an 
initial perturbation cannot stay coherent in a longitudinal 
spatial range spanning many wavelengths and hence many 
damping lengths. Instead, the system will select the frequen¬ 
cies and damping rates corresponding to the spatial periods 
enforced externally. 


Let us now focus on Fig. Eh, where one can see that for both 
density ratios, once in the leaky regime, k\ increases with de¬ 
creasing to and somehow levels off when to < 1 .5v^/a. On the 
other hand, k R decreases when to decreases from to c before in¬ 
creasing monotonically when to further decreases. Let to m denote 
the point where k R reaches a local minimum. It then follows that 
the apparent group speed dto/dk R < 0 for to < to m . One may 
question whether the sausage waves in this low-frequenc y por- 
tion can impart energy upwards. As discussed in Brillouinl (119601 
chapter V), when waves are heavily damped, the apparent group 


velocity may not represent the velocity at which energy propa¬ 
gates. In this case, we may directly use Eq. ([9]i to evaluate the 
z-component of the wave energy flux density. In view of Eq. ([5]i. 
one finds that 

(F z ) = |tv| 2 (—)exp(-2k : z), 

o7T \ tO ) 


which is positive for positive k R /w, meaning that the low- 
frequency waves in question can still direct their energy up¬ 
wards. However, in this case the plasma cylinder is such an inef¬ 
ficient waveguide that the wave energy is attenuated over a lon¬ 
gitudinal distance of the order of the cylinder radius. 

More insights can be gained by further comparing the disper¬ 
sive behavior of leaky standing and propagating sausage waves 
as given in Fig. E] To start, let us note that for standing waves 
with real k and complex to (propagating waves with real to and 
complex k) the radial energy flux density in the external medium, 
when averaged over a longitudinal wavelength 2n/k (a wave pe¬ 
riod 2k/ to), evaluates to 


27T i 

(Fr) — — J |/5tot,e 


Re 


itO P tot,e 
Pe Ptot,e 


,~2k[Z 


standing 

propagating. 


( 12 ) 


W e note that (F r ) for standing waves was originally derived 
in lCaTi\¥T986l. Eq. (3.2)). (See Appendix A for a derivation 
of Eq. CE3 .) Be the waves standing or propagating, p to t, e 00 
//‘j 1 \p c r) at large distances can be approximated by 



(13) 


resulting in p[ ote /ptot,e ~ ipe- It then follows that 


. 27T I |2 „ I to \ 

{ r) " If l /, t°t’ c l Re (^j 



standing 

propagating. 


(14) 


From Fig. Eb one can find that Im (p e ) < 0 for leaky stand¬ 
ing modes, meaning that |p tot>e | tends to grow exponentially 
(Eq. d. barring the r -1 ^-dependence) and an eigen-mode 
analysis does not allow us to examine the energetics of the sys¬ 
tem. However, this does not mean that this analysis is physi¬ 
cally irrelevant because what matters is that the apparent tem¬ 
poral damping can be accounted for by the outwardly going 
energy flux density (( F r ) > 0 because R e.(to/p e ) > 0, see 
Eq. d). Indeed, numerical studies starting with initial standing 
waves suggest that the temporal damping after a transient stage 
matches exac tly the attenuation rate given by the eigen-mode 
analysis dTerradas et al.ll2007l) . For propagating waves, one finds 
from Fig. Eh that once again Im (p j < 0 and R e(to/p e ) > 0. 
Hence similar to the standing case, an eigen-mode analysis does 
not permit an investigation into the energetics of propagating 
waves. However, if a coronal structure is perturbed with a har¬ 
monic boundary driver, one expects to see that the apparent spa¬ 
tial damping after some transient phase will be given by the 
attenuation length 1/ki obtained from this eigen-mode analy¬ 
sis. And this attenuation is once again associated with the out¬ 
wardly directed (F r ). We note that this expectation can be read¬ 
ily numerically tested, in much the same way that the expected 
spatial d ampi ng of kink waves due to resonant absorption was 
tested (iPascoe et al.ll201.3~l) . 
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4. Summary 

This study is motivated by the apparent lack of a dedicated study 
on the role that wave leakage plays in spatially attenuating prop¬ 
agating sausage waves supported by density-enhanced cylinders 
in the corona. To this end, we worked in the framework of cold 
magnetohydrodynamics (MHD), and numerically solved the dis¬ 
persion relation (DR, eq. (0)) for complex-valued longitudinal 
wavenumbers k - k R + ik t at given real angular frequencies to. 
To validate our numerical results, we also provided the analytical 
approximations to the full DR in the low-frequency limit to —* 0 
and in the neighborhood of to c , the critical angular frequency 
separating trapped from leaky waves. Our solutions indicate that 
while sausage waves can propagate for to < to c and can direct 
their energy upwards, they suffer substantial spatial attenuation. 
The attenuation length (1/ki) is of the order of the cylinder ra¬ 
dius a and shows little dependence on frequency or the density 
contrast between a coronal structure and its surroundings for 
to < 1.5i>Ai/a, where v\\ is the Alfven speed in the cylinder. This 
means that when a coronal cylinder is subject to boundary per¬ 
turbations with a broadband spectrum (e.g., granular motions), 
wave leakage removes the low-frequency components rather ef¬ 
ficiently. A comparison with the solutions to the DR for stand¬ 
ing waves (real k, complex to) indicates that a close relation¬ 
ship between propagating and standing waves exists only when 
the waves are trapped or weakly damped. In addition, while an 
eigen-mode analysis does not allow an investigation into the en¬ 
ergetics of propagating waves, the attenuation length is expected 
to play an essential role in numerical simulations where coronal 
structures are perturbed by harmonic boundary drivers. 
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Appendix A: A derivation of the averaged energy 
flux densities 

This section offers a derivation of the averaged energy flux densi¬ 
ties given in Eqs. ([8]), © and (fl2l) . following a procedure similar 
to that adopted by Li & Li : ( 20071 page 1225). Let us first con¬ 
sider standing waves for which the axial wavenumber k is real, 
but the angular frequency to is allowed to be complex-valued 
(to = ojr + ito\). Evaluating a perturbation 6f(r, z; t) with Eq. O 
at given values of [r, f] yields that 

6f( r, z; t) = Re {/ exp(ife)}, (A1) 

where 


/ = f(r) exp (-itot) = f(r) exp ( -ito R t ) e" ,f . (A2) 

With another perturbation 5g(r,z',t) in the same form, one finds 
that the product 6f6g averaged over a wavelength A - 2n/k is 


(6f6g) ( r , t) 

i r A 

- J 6f(r,r,t)5g(r,z\t)Az 

1 

- I Re {/exp(fc)} Re jg exp(/Lz)} dz 
1 W 

- I [/ R cos (kz) - fi sin(Lz)] [g R cos (kz) - gi sin(fo)] dz 

- J [/r0r cos 2 (fc) + figi sin 2 (kz) 

- (/r0i + Mr) cos (kz) sin(Lz)] dz 
^ (/r0r + Ml) 

xRe {fy*) - xRe (f*g) i (A3) 


where we have used the shorthand notations / R = Re/ and 
fi = Im/. By noting that / is expressible in terms of / through 
Eq. (IA2b . one finds that 

(SfSg)(r,t) = 1r e[f(r)g*(r)\e 2 ^ = ^Re \f(r)g(r)\ e 2 "Z(A4) 
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Now the energy flux density averaged over a wavelength can be 
evaluated with its definition, Eq. ([6]), the results being 


(F r ) = ^Re (ptotK) e 2 " lf - ^Re (Pi ot v r ) e 2 " 1 ', (A5) 

(A6) 

Furthermore, by noting that Eq. © allows v r to be expressed by 


v r 


. U ^Plot 

V B 2 ’ 


(A7) 


one finds that 


(Fr) 


2 71 

2 71 

^ Re 


. CO 

p tot\- l J£ p tot 




• ^ Aot 

AotAot I -l 2 ^ 

A* Ptot 


^2wr/ 


2?r IptotPn 

-^ Re 


\ B PtotJ 


(A8) 


This expression is valid both in and outside the cylinder. When 
applied to the external medium, it results in the first expression 
in Eq. (fl2l >. 

Now consider propagating waves for which io is real, 
whereas k is allowed to be complex-valued (k - k R + ik{). In 
this case evaluating a perturbation 6f{r, z; 1) with Eq. <□} at given 
values of [r, z] yields that 

<5/0, z; 0 = Re {/ exp(-/wr)}, (A9) 

where 


/ = /0) exp (, ikz ) = f(r) exp (ik R z) e hz . (A10) 


To evaluate the product 6f6g averaged over a wave period T = 
2i t/u>, one can follow the same procedure as in Eq. (IA3l i by re¬ 
placing kz with -cot, the result being 


(SfSg) (r, z) 

i r T 

= - 6f(r,z\t)6g(r,z;t)dt 

1 Jo 

= Irc (/§*)= ^Re(fg). (All) 

With the aid of Eq. (lAlOb which relates / to /, one finds that 

(6f6g) (r, z) = ^Re [/(r)sfO)] eT 2k ' z 

= ^Re \f( r )g(r)\e- 2k ' z . (A12) 

The energy flux density averaged over a period follows from the 
definition, Eq. ©, and the results are given by Eqs. © and ©. 
The second expression in Eq. (fT2l i. appropriate for propagating 
waves, can be derived in view of Eq. (fA7l 
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